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Abstract 

We have made a calculation of the propagation of positrons from dark-matter particle annihilation in the 

Galactic halo for different models of the dark matter halo distribution using our 3D code. We show that 

the Green's functions are not very sensitive to the dark matter distribution for the same local dark matter 

Qs I energy density. We compare our predictions with computed cosmic ray positron spectra ("background") for 

the "conventional" cosmic-ray (CR) nucleon spectrum which matches the local measurements, and a modified 

spectrum which respects the limits imposed by measurements of diffuse Galactic 7-rays, antiprotons, and 

f-H ! positrons. We conclude that significant detection of a dark matter signal requires favourable conditions and 

^ I precise measurements unless the dark matter is clumpy which would produce a stronger signal. Although our 

■ conclusion qualitatively agrees with that of previous authors, it is based on a more realistic model of particle 

^ . propagation and thus reduces the scope for future speculations. Reliable background evaluation requires new 

I accurate positron measurements and further developments in modelling production and propagation of cosmic 

^—* • ray species in the Galaxy. 

cn '. 1 Introduction: 

^P| I Investigations of galaxy rotation, big-bang nucleosynthesis, and large-scale structure formation imply that 

(^ ' a significant amount of the mass of the universe consists of non-luminous dark matter (Trimble 1989). Among 

On . the favored particle dark matter candidates are so-called weakly interacting massive particles (WIMPs), whose 

^ ! existence follows from supersymmetric models (see Jungman, Kamionkowski, & Griest 1996 for a review). 

i-^ I A pair of stable WIMPs can annihilate into known particles and antiparticles and it may be possible to detect 

I ■ WIMPs in the Galactic halo by the products of their annihilations. Though the microphysics is quite well un- 

P , derstood and many groups make sophisticated calculations of the spectra of annihilation products for numerous 

"^ I WIMP candidates which include many decay chains (e.g., Baltz & Edsjo 1998), there are still uncertainties 

C^ ■ in the macrophysics which could change the estimated fluxes of WIMP annihilation products by 1-2 orders 

^ . of magnitude, making predictions for their detection difficult. The most promising is perhaps the positron 

k>( I signal since it can appear at high energies where the solar modulation is negligible, but its strength depends 

V^ ■ on many details of propagation in the Galaxy. The "leaky box" model is often used (e.g., Kamionkowski & 

Turner 1991), a simplified approach which may not be applicable in the case of positrons. On the other hand, 

progress in CR positron measurements is anticipated since several missions operating or under construction 

are capable of measuring positron fluxes up to 100 GeV (e.g. experiments gas-RICH/CAPRICE: Barbiellini et 

al. 1997, and PAMELA: Adriani et al. 1997). Therefore, more accurate calculation of the positron propagation 

is desirable. 

We have developed a numerical method and corresponding computer code (GALPROP) for the calculation 
of Galactic CR propagation in 3D (for an overview of our approach and results see Strong & Moskalenko 
1999, and also papers OG 2.4.03, OG 3.2.18 in proceedings of this conference). Briefly, the idea is to develop 
a model which simultaneously reproduces observational data of many kinds related to cosmic-ray origin and 
propagation: directly via measurements of nuclei, antiprotons, electrons, and positrons, indirectly via 7-rays 
and synchrotron radiation. Here we use our model for calculation of positron propagation in different models 
of the dark matter halo distribution (Moskalenko & Strong 1999). To be specific we will discuss neutralino 
dark matter, although our results can be easily adopted for any other particle dark matter candidate. 
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2 Green's functions: 

The positron flux at the solar position is given by 
dF 



^ = Jde G{E, e) ^ BJ,{e) [cm-^ s'^ st~' GeV"!], 



(1) 



where /(e) is the source function which describes the spectrum of positrons from neutralino annihilation, 
G{E, e) is the Green's function for positron propagation in the Galaxy, and Bi is the branching ratio into a 
given final state i. The Green's function thus includes all details of the dark matter mass distribution, neutralino 
annihilation cross section, and Galactic structure (diffusion coefficient, spatially and energy dependent energy 
losses etc.). We can write it in the form: 



pI 



G{E, e) = {av) -^ g{E, e) [cm-^ g-^ sr-^ GeV" 
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where {av) is the thermally averaged annihilation cross section, po is the local dark matter mass density, m^ is 
the neutralino mass, and we have introduced a function g{E, e) which describes the positron propagation for a 
given dark matter mass density distribution 
in the halo. ^"2 g^ 

Following Kamionkowski and Kinkhab- 
wala (1998) we consider three different dark 
matter mass density profiles which match 
the Galactic rotation curve. The canoni- 
cal "isothermal" sphere profile, the spheri- 
cal Evans model, and an alternative model. 
For each given model we calculate the func- 
tion g{E, e) defined in Eq. Q, which gives 
the positron flux at the solar position corre- 
sponding to the positron source function in 
the form of a Dirac ^-function in energy. The 
positron propagation is calculated in a model 
which was tuned to match many available 
astrophysical data (Strong & Moskalenko 
1998, Strong, Moskalenko, & Reimer 1999). 
Since the halo size in the range Zh = A — 10 
kpc is favored by our analyses of B/C and 
^'^Be/^Be ratios and diffuse Galactic 7-ray 
emission, we consider two cases z/^ = 4 
and 10 kpc which provide us with an idea 
of the possible limits. The preferred neu- 
tralino mass range following from accelera- 
tor and astrophysical constraints is 50 GeV 
<m^< 600 GeV (Elhs 1998), and we con- 
sider positron energies e < 824 GeV which 
cover this range. 

Fig. [2 shows our calculated ^(-functions 
for different models of the dark matter dis- 
tribution: "isothermal", Evans, and alterna- 
tive. The curves are shown for two halo sizes 
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Figure 1 : Calculated (7-functions for different models of the 
dark matter distribution: (a) "isothermal", (b) Evans, (c) alter- 
native. Upper curves Zh = 10 kpc, lower curves Zh = A kpc, 
e = 1.03, 2.06, 5.15, 10.3, 25.8, 51.5, 103.0, 206.1, 412.1, 
824.3 GeV. The units of the abscissa are 10^^ GeV cm sr~^ 



Zh = 4: and lOkpc and several energies e = 1.03, 2.06, 5.15, 10.3, 25.8, 51.5, 103.0, 206.1, 412.1, 824.3 GeV. 
At high energies, increasing positron energy losses due to the inverse Compton scattering compete with the 
increasing diffusion coefficient, while at low energies increasing energy losses due to the Coulomb scattering 
and ionization (Strong & Moskalenko 1998) compete with energy gain due to reacceleration. The first effect 
leads to a smaller sensivity to the halo size at high energies. The second one becomes visible below ~ 5 GeV 
and is responsible for the appearance of accelerated particles with E > e. 

It is interesting to note that for a given initial positron energy all three dark matter distributions provide 
very similar values for the maximum of the (/-function (on the E'^g{E, e) scale), while their low-energy tails 
are different. This is a natural consequence of the large positron energy losses. Positrons contributing to the 
maximum of the (j(-function originate in the solar neighbourhood, where all models give the same dark matter 
mass density. The central mass density in these models is very different, and therefore the shape of the tail is 
also different since it is produced by positrons originating in distant regions. As compared to the isothermal 
model, the Evans model produces sharper tails, while the alternative model gives more positrons in the low- 
energy tail. At intermediate energies (-^ 10 GeV) where the energy losses are minimal, the difference between 
Zh = 4: and 10 kpc is maximal. Also at these energies positrons from dark matter particle annihilations in the 
Galactic center can contribute to the predicted flux. This is clearly seen in the case of the alternative model 
with its very large central mass density (Fig.fO;, Zh = 10 kpc). 

3 Positron fluxes: 

An important issue in the interpretation of the positron measurements is the evaluation of the "background", 
positrons arising from CR particle interactions with interstellar matter. Though the parameters of the propaga- 
tion and the Galactic halo size can be fixed in a self-consistent way using CR isotope ratios, the ambient CR 
proton spectrum on the Galactic scale remains quite uncertain. The only possibility to trace the spectrum of 
nucleons on a large scale is to observe secondary products such as diffuse 7-rays, positrons, and antiprotons. 

In order to show the effect of varying of the ambient proton spectrum, we compare our results with two 
models for the CR positron "background". These are a "conventional" model (model C) which reproduces the 
local directly measured proton and Helium spectra above 10 GeV (where solar modulation is small), and a 
model with modified nucleon spectrum (model HEMN), 
which is flatter below 20 GeV and steeper above, and 
results from our analysis of Galactic diffuse 7-ray emis- 
sion. The "background" spectra are slightly dependent 
on the halo size. Since all secondary particles are pro- 
duced in the Galactic plane, increasing the halo size re- 
sults only in a small decrease of the flux at high energies 
due to larger energy losses. The propagation parame- 
ters for these models are given in Strong & Moskalenko 
(1998) and Strong, Moskalenko, & Reimer (1999), and 
the formalism for calculation of secondary positrons is 
described in Moskalenko & Strong (1998). 

We do not intend to make sophisticated calcula- 
tions of positron spectra resulting from numerous de- 
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Figure 2: Our predictions for two CR positron 
"background" models (C and HEMN: heavy solid 
cay chains such as best done by, e.g., Baltz & Edsjo lines), and positron signals from neutralino annihila- 
(1998) for many WIMP candidates. Instead, for illustra- tion for m^ = 5.15, 10.3, 25.8, 103.0, 206.1, 412.1 
tion purposes, we simplify our analysis by treating the (thin solid lines): (a) z^ = A kpc, (b) z^ = 10 
annihilation to W^ and Z°-pairs. For ruy. < niw we kpc- In the case of m^ = 103.0 GeV, the signal 
consider only the direct annihilation to e+e" pairs. In Plus background (model C) is shown by the dotted 
the first case we use the cross sections for a pure Hig- line. Data and the best fit to the data (dashes) are 
gsino (Kamionkowski & Turner 1991), in the latter case from Barwick et al. (1998, HEAT collaboration). 



we take B • (av) = 3 x 10"^^ cm^ s~^ and monoenergetic positrons. These parameters can be considered as 
optimistic, but possible. To maximize the signal we further choose the Galactic halo size as 10 kpc. 

Fig. |2l shows our predictions for the two CR positron "background" models together with HEAT data 
(Barwick et al. 1998) and positrons from neutralino annihilation. It is seen that the predicted signal/background 
ratio has a maximum near m^ ~ mw, while even in the "conventional" model the background is nearly 
equal to the signal at its maximum. It is however interesting to note that our calculations in this model show 
some excess in low energy (< 10 GeV) positrons where the measurements are rather precise but the solar 
modulation is also essential. If this excess testifies to a corresponding excess in interstellar space and if the 
positron background correspond to our "conventional" calculations, it could be a hint for the presence of the 
dark matter (Baltz & Edsjo 1998, Coutu et al. 1999). Our HEMN model fits the HEAT data better (no excess) 
and thus provides more background positrons. (This shows that in principle a good fit to positron data, which 
is consistent also with other measurements such as 7-rays and antiprotons is possible without any additional 
positron source.) Under such circumstances a significant detection of a weak signal would require favourable 
conditions and precise measurements. Though this conclusion qualitatively agrees with that of Baltz and Edsjo 
(1998) and several earlier papers, it is based on a more realistic model of particle propagation and thus reduces 
the scope for future speculations. 

4 Conclusions: 

Our propagation model has been used to study several areas of high energy astrophysics. We use this 
model for the calculation of positron propagation in different models of the dark matter halo distribution. We 
have shown that the Green's functions are not very sensitive to the dark matter distribution for the same local 
dark matter energy density. This is a natural consequence of the large positron energy losses. We compare 
our predictions with the computed CR positron "background" for two models of the CR nucleon spectrum. A 
correct interpretation of positron measurements requires reliable background calculations and thus emphasizes 
the necessity for further developments in modelling production and propagation of CR species in the Galaxy. 
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